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ABSTRACT 

The strength and spectrum of the extragalactic magnetic field are still unknown. 
Its measurement would help answer the question of whether galactic fields are purely 
a primordial relic or were dynamically enhanced from a much smaller cosmological 
seed field. In this letter, we show that the composition, spectrum, and directional 
distribution of extragalactic ultrahigh energy cosmic rays with energies above ~ 10 18 eV 
can probe the large scale component of the extragalactic magnetic field below the present 
observational upper limit of 10 -9 Gauss. Cosmic ray detectors under construction or 
currently in the proposal stage should be able to test the existence of the extragalactic 
magnetic fields on scales of a few to tens of Mpc and strengths in the range ~ 10 -10 — 
10~ 9 Gauss. 

Subject headings: cosmic rays — magnetic fields — gamma rays: theory 

1. Introduction 

Magnetic fields are universally present in astronomical bodies ranging from the Earth to the 
distant quasars, but it is still unknown if magnetic fields permeate the universe as a whole. Astro- 
physical magnetic fields may arise due to the existence of a primordial magnetic field that grows 
as galaxies form. The discovery of an extragalactic magnetic field on scales larger than virialized 
systems (i.e., larger than clusters of galaxies) would reveal the presence of a primordial field. The 
existence of such a primordial field would help understand the origin of astrophysical magnetic 
fields and may open a new window into processes ocurring in the early universe. 
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In this letter, we show that the study of ultra high energy cosmic rays (UHE CRs) with energies 
above ~ 10 18 eV can probe primordial fields below the current upper bound. Our proposed method 
is complementary to traditional ones (e.g., Kronberg 1994) and more recent suggestions (e.g., 
Plaga 1995). At present, the most stringent constraint on large scale extragalactic fields comes 
from limits on the Faraday rotation of light coming to us from distant quasars. The upper bound 
on a widespread, all-pervading field is ~ 1CT 9 Gauss (e.g., Kronberg 1994). There are weaker 
constraints derived from the synchrotron emission from nearby galaxy clusters (Kim et al. 1989) 
and the cosmic microwave background isotropy. 

UHE CR nucleons from extragalactic sources are attenuated in energy while propagating 
through the cosmic microwave background (CMB). Above a few 10 19 eV, nucleons produce pi- 
ons on the CMB photons and the energy of the cosmic ray nucleons is degraded rapidly, which is 
known as the Greisen-Zatsepin-Kuz'min (GZK) effect (Greisen 1966; Zatsepin & Kuz'min 1966). 
Produced pions, on the other hand, decay into secondary leptons such as electrons, muons, and 
neutrinos, and photons. Since muons decay to electrons and neutrinos, the final secondary particles 
are photons, electrons, and neutrinos, among which photons are more readily detectable. 

Very energetic secondary photons and electrons couple to form an electromagnetic (EM) cas- 
cade. The 7-rays produce electron pairs on the CMB and the universal radio background photons 
(see Lee & Sigl 1995, Lee 1995 and references therein for more detailed discussions). The result- 
ing electrons (or positrons) in turn upscatter background photons via inverse Compton scattering 
(ICS), thus completing a cycle. It is through these two processes that an EM cascade develops 
in the intergalactic medium. As a result, if one propagates a purely protonic spectrum, one gets 
a processed nucleon spectrum with the GZK cutoff, and a secondary EM (photons and electrons) 
spectrum. 

The extragalactic magnetic field (EGMF) influences the UHE CR flux mainly through their 
charged components, namely the primary hadrons and the secondary electrons. In the energy 
range under consideration the hadrons are deflected with negligible energy loss due to synchrotron 
radiation whereas the electrons are negligibly deflected before they lose most of their energy. In §2, 
we explore the effect of the electron energy loss due to the EGMF on the secondary EM cascade 
spectrum. We then discuss the deflection of the hadronic component in the EGMF, in §3. Finally 
in §4, we summarize our findings. 



The EGMF plays a crucial role in the development of the EM cascade. In the presence of the 
magnetic field, electrons lose energy via synchrotron radiation loss, which is given by 



2. Extragalactic Magnetic Field and Ultrahigh Energy 7-Ray Flux 
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where B is the strength of the large scale EGMF, r$ is the classical electron radius, and m e is the 
electron mass. 

In fig. 1, we show the rates of ICS and synchrotron loss for electrons. Whereas the rate of ICS 
responsible for EM cascade development decreases with energy, the synchrotron loss rate increases 
with energy. Therefore, in a narrow energy range a transition occurs between a regime where ICS 
is dominant and electrons couple to photons efficiently and another where electrons are rapidly lost 
due to synchrotron loss and the cascade is suppressed. Below ~ 10 20 eV (the threshold for pair 
production on the radio background), cascade development in the absence of an EGMF would give 
rise to a generic power law photon spectrum with index ~ —1.5. The above mentioned transition 
in the secondary 7-ray spectrum will therefore occur between this generic cascade shape and a 
synchrotron loss dominated spectrum. As long as synchrotron quanta can be neglected the latter is 
given by the photons produced "directly" by source injection or from pion production by nucleons, 
before undergoing pair production in the low energy photon background. 

In a magnetic field of strength B measured in G (Gauss) the synchrotron spectrum produced 
by an electron of energy E e peaks at 

E - = 6 - 8xlol3 (loitv) 2 (lFb) eV - (2) 

and falls off exponentially at higher energies. In the following we assume that the observable EM 
flux is energetically dominated by 7-rays and electrons with energy E <^ 10 21 eV. Then, according 
to eq. |2]|, the contribution of synchrotron radiation to the 7-ray flux above 10 18 eV can be safely 
neglected as long as B <^ 10 -6 G everywhere. 

The energy where the transition in the 7-ray spectrum occurs is in general a function of the 
magnetic field strength and the background photon spectrum, but not a function of the source 
distance or the injection spectrum. If we consider only the CMB, the relation between the tran- 
sition energy E tr and the magnetic field strength is Et T oc B . In order to include the less well 
known diffuse extragalactic radio background into consideration we adopt its usual description by 
a power law with an overall amplitude and a lower frequency cutoff as parameters (Clark, Brown, 
&i Alexander 1970) the latter one being the main source of uncertainty. Using a cutoff at 2 MHz 
as suggested by Clark et al. (1970) the above relation is modified to 

Eu * 10 19 (j^g) 13 eV a 10_10 G ) • ( 3 ) 

For the same magnetic field a cutoff at lower frequencies would increase the rate of ICS of the then 
more abundant low frequency radio photons and thus the value for Et r (see fig. 1). Assuming that 
the radio cutoff frequency lies somewhere in the range between 0.5 MHz and 3 MHz, for a given 
Et T the EGMF strength B is uncertain within about a factor 5. 

Therefore, for B ^ 10~ 10 G it is possible to approximately determine the EGMF strength by 
searching for a dip in the 7-ray flux below 10 20 eV which would mark a transition between an ICS 
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and a synchrotron loss dominated regime. For B <^ 10 G, this transition occurs above the pair 
production threshold on the radio background where the 7-ray flux increasingly depends on several 
unknown factors such as the charged cosmic ray flux above the GZK cutoff. Thus, even though the 
7-ray flux can be comparable to the nucleon flux above 10 20 eV, a discussion of possible magnetic 
field signatures in its spectrum would presently be too speculative. 

We developed a numerical code for the propagation of nucleons, photons, and electrons through 
the intergalactic medium which employs a transport equation formalism, the details of which can 
be found in Lee (1995). The observed UHE CR flux below 10 20 eV (see, e.g., Bird et al. 1994; 
Yoshida et al. 1995) is reproduced quite well by a diffuse distribution of sources injecting protons 
with a spectrum cx E~ 2 ' 3 up to some maximal energy considerably beyond the GZK cutoff (Yoshida 
& Teshima 1993; Sigl et al. 1995). For the calculations presented here we therefore adopted this 
proton injection spectrum with a maximal energy of 10 22 eV (see figures). The EGMF enters the 
calculation via the synchrotron loss of electrons. 

The transition between ICS and synchrotron loss domination can be easily seen in fig. 2, which 
shows the processed nucleon and photon spectra for a single source at a distance of 30 Mpc for a 
range of EGMF strengths. In general, one expects a distribution of cosmic ray sources rather than 
a single source at a fixed distance. In fig. 3, we show the diffuse spectrum from a continuous source 
distribution extending up to 1 Gpc. We assume a flat universe with zero cosmological constant and 
a Hubble constant of Hq = 75 km sec _1 Mpc _1 , and a comoving source density scaling as (1 + z) 2 
in redshift z as in some "bright phase" models of CR sources (e.g., Yoshida & Teshima 1993; Hill 
& Schramm 1985 and references therein). The results are not very sensitive to these choices. In 
the diffuse case, the 7-ray to nucleon flux ratio tends to be smaller than for a single source, and 
the spectral features are not as pronounced, but still detectable for B ^ 10~ 10 G. 

The "extragalactic magnetic field" in this analysis refers to the average component of the 
EGMF normal to the line of propagation. Primordial magnetic fields are expected to have very 
little structure on scales below ~ few Mpc (Jedamzik, Katalinic, & Olinto 1995), but condensed 
structures such as galaxies and clusters of galaxies can "polute" the intergalactic medium with 
stronger magnetic fields on smaller scales. Fortunately, the effect of the EGMF on the 7-ray 
spectral shape discussed here is most sensitive to the average field on large scales. When the EM 
cascade goes through a strong field region, the electrons in the cascade lose energy rapidly and drop 
out of the UHE range. In contrast, the UHE nucleon and 7-ray fluxes are usually hardly affected 
directly by the radiation field of the intervening object (Stecker et al. 1991; Szabo & Protheroe 
1994; Norman, Melrose, & Achterberg 1995). After escaping the object, the cascade redevelops 
and the cascade spectrum recovers quickly at only a slightly smaller amplitude. The influence of 
intervening objects decreases with increasing strength of the large scale field and becomes important 
only when the objects are very close to the observer (e.g., less than ~ 5 Mpc away), or when their 
linear size is significant compared to the total propagation distance. We can derive conditions for 
the filling factors for such objects by requiring that they are much more sparsely populated along 
the line of sight to the source than the typical cascade regeneration length s c (e.g., ~ 5 Mpc). For 
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an average linear size k v of the intervening objects their filling factor f w must satisfy 

/iv « (4) 

s c 

For field galaxies the above relation is f g <C 10~ 3 , and for clusters of galaxies f c <C 0.1. The 
actual filling factors for galaxies, f g <^ 10~ 7 , and galaxy clusters, f c <^ 10~ 4 (Kolb & Turner 
1992; Nichol, Briel, & Henry 1994) satisfy these constraints. Therefore, intervening objects do not 
modify the above discussion substantially. One interesting exception may be the effect of nearby 
large structures such as the Virgo Cluster. If Virgo has strong magnetic fields, e.g. of order 10 -7 G 
on Mpc scales, the UHE 7-ray flux from background sources might be modified across Virgo's 
angular extension (see fig. 2). 

The ability of future detectors to study EGMF features crucially depends on the 7-ray to 
nucleon flux ratio. For nearby strong sources, the secondary 7-ray flux could be measurable with 
instruments which are sensitive to ratios down to ~ 1%. This could be achieved by the proposed 
Pierre Auger Project (see e.g., Boratav et al. 1992), which would also allow an angular resolution 
of ~ 1°. The case of a diffuse source distribution is more challenging; the 7-ray flux is typically 
smaller and the EGMF feature is less pronounced than for a single source at moderate distances 
(see figs. 2 & 3), but the dip in the 7 spectrum may still be detectable. Thus, measuring the 7-ray 
flux between ~ 10 18 eV and ~ 10 20 eV has the potential to either detect or find strong evidence 
against an EGMF £ 1(T 10 G. 



3. Charged Cosmic Ray Deflection by Extragalactic Magnetic Fields 

Here, we discuss the influence of the EGMF on the charged UHE CR flux from discrete sources. 
We restrict ourselves to the case of small deflection angles (for the opposite limit see, e.g. Wdowczyk 
& Wolfendale 1979; Berezinskii, Grigor'eva, & DogeP 1989). In this case, the energy spectrum of 
charged UHE CRs from a given source is not significantly altered as compared to a straight-line 
propagation. However, if the sources are strong enough to cause an anisotropy in the UHE CR 
flux, the directional correlation of "hot spots" with possible sources will depend on the EGMF. 
The following discussion relates to this anisotropic component of the charged UHE CR flux. 

As in §2, let us assume that the large scale EGMF can be characterized by a typical field 
strength B and a coherence length l c . Furthermore, we assume for the moment that the source 
distance r is smaller than the energy attenuation length A = E(dE/dr)~ 1 for a charged cosmic ray 
of energy E which can then be treated as approximately constant throughout propagation. For 
nucleons, A ~ 10 Mpc above the GZK cutoff (at E ~ 6 x 10 19 eV), and A ~ 1 Gpc much below 
the GZK cutoff. A more sophisticated analysis would require a Monte Carlo simulation of both 
UHE CR propagation and deflection. However, since data on both UHE CR and the EGMF are so 
sparse to date, we feel that a qualitative discussion of the principle effects is more appropriate at 
the moment. We now consider two cases. 
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(i) The source distance is smaller than the coherence length, r <^l c . Then, in vectorial notation, 
the deflection angle a is given by 

" - T — ^ (u£v) " Go**) (iPo) ( f * 6 ) • < 5 > 

where r is the radius vector pointing to the source, Z is the charge of the UHE CR component, 
f = r/|r|, and B = B/|B|. Thus, a correlation between the UHE CR flux of charge Z and energy 
E and source counterparts, systematically shifted by an angle a, would indicate that Z c <; r and for 
a known source distance r would allow to measure the combination B(r x B). The characteristic 
E- and Z dependence of a would provide an additional test for the hypothesis that the deflection 
is caused by an EGMF. 

(ii) The source distance is considerably larger than the coherence length, r>l c . In this case 
the deflection angle undergoes a diffusion process during propagation and the source shape in the 
UHE CR flux will be smeared out over a typical angle 

arms ^~ (WJ 1 / 2 = 1.1° Z (-^j (^-) ' (^) ' (^) . (6) 

Therefore, if sources appear spread out in the UHE CR flux of charge Z and energy E by a typical 
angle a, this would indicate that l c <; r and for a known source distance r would allow to measure 
the combination Bl\^ 2 . 

If the source distance is larger than the energy attenuation length, r ^ A, eqs. || and || tend 
to overestimate a. In fact, in the limit r S> A the deflection angle a "saturates" as a function of r 
and for approximately energy independent A, r has to be substituted by A and A/ 2 in eqs. || and 
||, respectively. 

Secondary 7-rays produced by the interactions of the charged UHE CRs are also expected to 
correlate with the sources. Due to their continuous production they will be smeared out over angles 
which are typically somewhat smaller than given in eqs. 

§ and f§. 

Sources which can act as suitable probes for the EGMF via the effects discussed above have 
to obey the following conditions apart from producing a detectable anisotropic UHE CR flux 
component: Their apparent angular size should be smaller than the deflection angle a. The same 
pertains to the apparent angular radius of a possible high magnetic field region around the source if 
it can cause deflections in excess of a. For example, a 1CP 6 G field over a scale ;> 100 kpc is possible 
in galaxy clusters (see, e.g., Kronberg 1994) and would completely bend around a 10 20 eV proton. 
However, as long as a detectable proton flux emerges from such an object and the above conditions 
are fulfilled, it could still be a suitable probe of the EGMF. Finally, there should be no intervening 
high magnetic field regions between source and observer which could cause bending by more than 
a. For example, for r ^ l c , A, this corresponds to the condition k v B- lv (\l c ) 1 ^ 2 B for linear scale 
and strength B n of the intervening field. This condition could well be satisfied along most lines 
of sight since known objects with high field regions like galaxies and galaxy clusters have a small 
filling factor / <; 10~ 5 . 
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In light of these conditions we believe that some of the nearby galaxy clusters and powerful 
field radio galaxies could well be suitable EGMF probes since they are expected to contribute 
significantly to the UHE CR flux (Rachen, Stanev, & Biermann 1993). The accuracy to which 
the EGMF bending can be determined is limited by the (to date) unknown additional bending 
by the galactic magnetic field of strength B g and scale height l g . Thus, according to eq. [||, the 
sensitivity of deflection measurements of the EGMF is restricted to field parameters satisfying 
Br ^ 10~ 9 G Mpc (Bg/fxG) (Z 9 /300pc) where the fudge factors are the parameter values usually 
assumed for the galactic magnetic field. 

4. Conclusions 

We discussed how composition, spectrum, and directional distribution of UHE CR above ~ 
10 18 eV can be used to gain information about the large scale (a few to tens of Mpc) EGMF. Spectral 
features in the 7-ray flux are sensitive to field strengths in the range ~ 10 -10 — 1CP 9 G. In a similar 
range, correlations between an anisotropic charged UHE CR flux component and possible sources 
could provide independent information on the EGMF including its polarization. Both effects should 
yield consistent estimates for the EGMF strength. Strong discrete sources detected in UHE CRs by 
future instruments with an angular resolution of 1° or better and a sensitivity to 7-ray to nucleon 
flux ratios of 1% or smaller would provide the best conditions for detecting an EGMF in the range 
~ 10 -10 — 10~ 9 G. Since these conditions are not unreasonable, UHE CRs have the potential to 
provide important information on properties and origin of the EGMF. 
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Fig. 1. — The rate of inverse Compton scattering (ICS; solid lines) and the fractional synchrotron 
loss rate (dashed lines) for different magnetic fields as functions of energy. 1 is for a frequency 
cutoff at 2 MHz in the diffuse radio background, and 2 is for 0.5 MHz, respectively. 
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Fig. 2. — Cosmic ray nucleoli spectrum (solid line) from a single source at a distance of 30 Mpc, and 
secondary 7-ray spectra for different magnetic fields. The short dashed lines are with no intervening 
objects, and the long dashed lines with dots are with a 5 Mpc sized object located at 12.5 Mpc 
from the observer. The magnetic field values are B = 0, 10~ 10 , 3 x 10~ 10 , and 10~ 9 G, respectively 
from the top. 
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Fig. 3. — Cosmic ray nucleon spectrum (solid line), and secondary 7-ray spectra (dashed lines) for 
a continuous source distribution up to lGpc distance whose comoving density scales as (1 + z) 2 . 
The sum of the nucleon and 7-ray flux is fitted to the combined data from the monocular Fly's Eye 
(Bird et al. 1994) and the AGASA experiment (Yoshida et al. 1995). Shown are la error bars, and 
84% confidence level upper limits. The magnetic field values are the same as in fig. 2. 



